Experimental Evidence of Time Delay Induced Death in Coupled Limit Cycle 

Oscillators 
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Experimental observations of time delay induced amplitude death in a pair of coupled nonlinear 
electronic circuits that are individually capable of exhibiting limit cycle oscillations are described. In 
particular, the existence of multiply connected death islands in the parameter space of the coupling 
strength and the time delay parameter for coupled identical oscillators is established. The existence 
of such regions was predicted earlier on theoretical grounds in [Phys. Rev. Lett. 80, 5109 (1998); 
Physica 129D, 15 (1999)]. The experiments also reveal the occurrence of multiple frequency states, 
frequency suppression of oscillations with increased time delay and the onset of both in-phase and 
anti-phase collective oscillations. 

PACS numbers: 05.45.Xt, 87.10. +c 



Coupled limit cycle oscillator models have been exten- 
sively studied in recent years because of the useful in- 
sights they provide into the collective behaviour of many 
physical, chemical and biological systems [0-|4). One of 
the simplest such models, the so-called Kuramoto model 
|E| , which retains only the phase information of each os- 
cillator and is valid in the limit of weak mutual cou- 
plings, displays a spontaneous transition to a synchro- 
nized collective state of a single frequency when the cou- 
pling strength exceeds a critical value. Similar collective 
behaviour is observed in many natural systems, such as 
the synchronous flashing of fireflies, the phase locking of 
cardiac pacemaker cells, and the collective chirping of 
crickets (i). Phase locking has also been demonstrated 
experimentally in arrays of coupled nonlinear electronic 
circuits ||. When the coupling becomes stronger, am- 
plitude effects become important and give rise to other 
interesting collective states, such as that of amplitude 
death in which the various oscillators pull each other off 
their periodic states and collapse to a state of zero am- 
plitude. The condition for such a state to occur is for 
the oscillators to have a broad dispersion in their natu- 
ral frequencies and for the coupling strength to exceed 
a threshold value. Thus, as has been pointed out in a 
number of theoretical studies [||, a collection of identi- 
cal limit cycle oscillators cannot display amplitude death. 
However, more recent investigations [Q indicate that the 
presence of finite propagation time delays in the coupling 
removes this restriction and predict the possibility of in- 
ducing the death state even in a system of two coupled 
identical limit cycle oscillators. Time delay is ubiquitous 
in most physical systems due to finite propagation speeds 
of signals, finite chemical reaction times, finite response 
times of synapses, etc., and its influence on the collective 
dynamics of coupled systems can have wide-ranging im- 
plications ||. It is important therefore to establish the 
experimental feasibility of such a death phenomenon. 



In this Letter we present experimental observations on 
time delay induced death in two coupled nonlinear cir- 
cuits that are individually capable of exhibiting limit cy- 
cle oscillations. A specially designed digital delay line 
gives precise control over the delay time in the coupling 
and permits us to explore a large area of parameter space. 
Observations on the phenomenon of amplitude death have 
recently been reported Q for a pair of optothermal oscil- 
lators that are thermally coupled and for which the oc- 
currence of death for strong couplings and its relation to 
Hopf bifurcations of the uncoupled and coupled systems 
have been experimentally verified. However, as the au- 
thors themselves point out, their experimental results are 
not conclusive about the role of delay in the death phe- 
nomenon and in particular they have not investigated the 
phenomenon of death islands or, indeed, that of the mul- 
tiplicity of death islands predicted by the theory of time 
delay induced death jjj. In our experiments we provide 
clear evidence of time delay induced death islands and 
their multiple connectedness in the parameter space de- 
fined by the coupling strength, time delay and frequency, 
for the case of two coupled identical oscillators. We also 
find the existence of multiple frequency states, frequency 
suppression of oscillations with increased time delay and 
the onset of both in-phase and anti-phase collective os- 
cillations. 

Our experimental system, schematically shown in 
Fig. 0, consists of two nonlinear LCR circuits coupled 
through a digital delay line (DDL). The individual oscil- 
lator circuits are a variant of the so-called Chua circuit, 
which has been widely employed in a number of nonlinear 
dynamical studies |lC| ]. The nonlinear resistive element 
i?jv has the typical V — I characteristic, with negative re- 
sistance around the origin and positive slopes away from 
the origin as shown in Fig. |^. This enables each indi- 
vidual oscillator circuit in the uncoupled state to sustain 
limit cycle oscillations with a characteristic frequency of 
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1 / y/LC. The two channel DDL taps the input signal at a 
certain rate using an analog to digital converter (in this 
case an ADC0809 with a conversion rate of 100/is) and 
stores it in a random access memory (RAM) bank. The 
stored bits are read using a logical circuit and converted 
back to an analog signal to be fed back into the circuit. 
The coupling between the oscillators is linear, resistive 
and proportional to the difference in the signal strengths 
of the two oscillators with a time delay. The coupling 
strength is varied by changing the resistances (Rk) that 
join the two oscillators. By using the current summation 
rules at nodes V\ and V2, a theoretical description of the 
circuit can be given as, 

Vi + g(Vi) Vi + ujI Vi = Ki[V s (t - r) - V(t)], (1) 

where i,j = l,2,i ^ j and g(Vi) = (l/C;)(-0.95 + 
0.66^ 2 + 322.0V" l 4 ) x 1CT 5 Amp/F, the nonlinear damping 
factor, is obtained from an empirical fit of the V — I char- 
acteristics of the nonlinear resistive element as shown in 
Fig. ||. The frequencies are u>j = 1/ yjLj Cj and the cou- 
pling strengths Kj = 1/ (CjRx )■ For our experiments 
on identical oscillators we have fixed C\ = C2 = 0.1/iF 
and varied L = L\ = Li and R = Rx t — Rk 2 such that 
u) = oji = L02 spans the range of 100s -1 to 1000s -1 . 




FIG. 1. The circuit diagram of two delay coupled limit cy- 
cle oscillators consisting of a C, L and a nonlinear resistance 
Rn connected in parallel. The OP-AMPs are buffer ampli- 
fiers and ideally draw no input currents. The digital delay 
line outputs V{ 2 = Vi,i(t — t). 
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FIG. 2. The V-I characteristics of the nonlinear component 
Rn- The continuous line is a polynomial fit of the experimen- 
tal points. 



The coupling strength varies in the range of 10s _1 to 
10 6 s -1 . Examples of our experimental results for the 
time evolution of the oscillator voltages as a function of 
the delay parameter are shown in Fig. || for K = 1000s -1 
and to — 837s -1 . The identical oscillators acquire an in- 
phase locked state (for r = 0.514 ms), death (for r = 2 
ms) and an anti-phase locked state (for r = 4.428 ms). 
A series of such observations are recorded for a fixed fre- 
quency and at different fixed values of K by varying the 
time delay parameter. The output voltage V of the os- 
cillator is monitored on the oscilloscope up to the point 
where it shows a sudden drop in value by several orders of 
magnitude, i.e., nearly to zero. The value of r is further 
increased up to the point where the signal from the oscil- 
lator again suddenly revives to a large value. An entire 
parameter space in K — t for amplitude death is thereby 
explored. In Fig. ^ the solid squares are the experimental 
data points marking the loci of the death regions. They 
clearly reveal closed death island regions in the plane of 
the coupling strength, K , and the time delay, t, for the 
common frequency of u> = 837s -1 . 




FIG. 3. Oscilloscope traces of the temporal behavior of the 
voltages of the two delay coupled oscillators for different delay 
times with K = 1000s -1 and lo — 837s -1 . Time is marked 
in units of 5 ms along the horizontal axis. The oscillators 
show (top panel) in-phase locking for r = 0.514 ms, (middle 
panel) amplitude death for r = 2.0 ms, and (bottom panel) 
anti-phase locking for r = 4.428 ms. 
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FIG. 4. Multiply connected death islands in K — r space 
for lu = 837s" 1 . Solid squares mark the experimental points 
while the solid and dashed lines are theoretical curves ob- 
tained from model equation (1). 

The solid and dashed lines are analytic curves obtained 
from a linear stability analysis of the origin for the model 
equation (Q). These curves are derived in a standard 
manner H from the characteristic equation for the eigen- 
value A, namely A 2 (A - % + K) 2 + 2w 2 (A - % + K)\ + 
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where A = J '^K - (^-) 2 , a t = 0.95 x 10 -5 and 
n = 1,2,.... The region enclosed by the intersection 
of Ti(l, K) and T2(l, K) is the primary death island and 
that enclosed between t\(2,K) and t%(2,K) is the sec- 
ondary (higher order) death island. We observe that the 
experimental observations not only display topological 
similarity to the analytic predictions but are also in good 
quantitative agreement with them. The deviations ob- 
served are primarily due to two reasons. First, the ana- 
lytic domains are sensitive to the values of the coefficients 
used in fitting the V — I characteristic of the nonlinear 
resistive element. The relative downward shift and larger 
width of the experimental islands compared to the ana- 
lytic region are due to these estimation errors. Second, 
the deviation at very high values of K (and hence low 
Rk) is related to limitations on the maximum current 
that can flow through Rk. for the rated biased voltages 
of our experimental circuit. The existence of the higher 
order island is particularly significant since it confirms 
the theoretical finding of multiple connectivity of the sta- 
bility region of the origin |t]] . 
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FIG. 5. Amplitude death region in the plane of (lu, t) for a 
fixed coupling strength of K — 500s -1 . 



In Fig. [jj, the experimental death region is shown at 
a constant coupling strength K = 500s -1 for different 
values of the frequency. The death region ceases to exist 
below a certain threshold frequency, in accordance with 
theoretical predictions @. On either side of this region, 
there are phase locked states exhibiting in-phase locked 
oscillations for small time delays and anti-phase locked 
oscillations for longer time delays. 

Another important characteristic of time delay sys- 
tems in general [[ll],[l2), and delay coupled limit cycle 
oscillators in particular is the existence of multi- 
ple phase locked states. When the oscillators are iden- 
tical, these phase locked states coincide with in-phase 
(phase difference, = 0) and anti-phase (cf> — ir) so- 
lutions. In Fig. ^(a), the domains of existence of these 
solutions in t space are indicated. Here /„ represent 
the in-phase solutions and A n represent the anti-phase 
solutions. These domains have been determined exper- 
imentally. For a fixed r the various existing multiple 
states can be accessed with slight perturbations to the 
system, i.e., by slightly changing the initial conditions. 
We have also observed hysteretic behaviour between the 
various branches. In Fig. ^(b), the frequencies of the first 
few multiple states are plotted for increasing values of r. 
Note that the collective states can be much higher than 
the individual oscillator frequencies and that 
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FIG. 6. (a) Coexistence of in-phase and anti-phase locked 
states, and (b) suppression of the phase locked states as r is 
increased for K = 1000s -1 and lu = 837s -1 . 
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their magnitude decreases with increasing values of time 
delay. Such a phenomenon of time delay induced fre- 
quency suppression has been discussed earlier in the 
context of phase only oscillators . 

In conclusion, we have experimentally demonstrated 
for the first time the existence of delay induced am- 
plitude death in a pair of coupled identical nonlinear 
oscillators R], We have also demonstrated the coex- 
istence of in-phase and anti-phase states with multiple 
frequencies and the suppression of the common frequency 
with time delay. The observation of multiply connected 
death islands (for identical oscillators) in the param- 
eter space of the coupling strength and time delay is 
our most important result. This result not only bears 
topological resemblance to earlier model calculations , 
but also shows good quantitative agreement with ana- 
lytic estimates obtained from a theoretical model of the 
experimental circuits. It should be noted that our exper- 
imental system of coupled nonlinear electronic circuits 
(and its theoretical model) is capable of richer nonlin- 
ear behaviour and goes beyond the simple normal form 
model used in previous theoretical studies 0. In this 
sense the results on delay induced amplitude death can 
be said to have a broader validity and therefore wider 
applicability to physically realistic systems. 

We now briefly discuss a few interesting directions for 
future experiments that are relevant to our present re- 
sults as well as to practical applications. A natural next 
step is to couple more than two oscillators to explore 
some of the other theoretical predictions discussed in |Q . 
The principal technical challenge and the major expen- 
sive elements for such an endeavor are the construction of 
additional delay lines with enhanced number of channels. 
Such an effort is presently underway. Another interesting 
possibility is to move away from the limit cycle regime 
and investigate the effect of time delay on the synchro- 
nization of coupled chaotic oscillators - an area of much 
current scientific interest Jl3| . In principle this should be 
possible with our present experimental system since the 
basic Chua circuits are capable of displaying chaotic oscil- 
lations. It would be necessary however to employ a faster 
delay line. Time delay could facilitate phase synchroniza- 
tion in view of its significant influence on phenomena like 
phase slips |Q. A final area of practical interest is the 
stability of these limit cycle collective states to external 
noise. Past numerical studies indicate that time delay 
induced higher frequency states can exhibit metastable 
behaviour and decay to the lowest periodic state in the 
presence of sufficient external noise [0|. We do not ob- 
serve this behaviour since by design our system has very 
low intrinsic (thermal) noise in order to facilitate access 
to the higher frequency states. It would be interesting 



to introduce a suitably designed noise source in the in- 
dividual oscillator circuits to experimentally test these 
theoretical predictions including the influence of noise on 
the threshold for onset of synchronization |l5| . 

We acknowledge the important contribution of H. S. 
Mazumdar in the design, development and construction 
of the digital delay line. 
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